Magnetism in small bimetallic Mn-Co clusters 
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Effects of alloying on the electronic and magnetic properties of Mn^Co^ (x + ^=n=2-5; x=0- 
n) and Mn2Coii nanoalloy clusters are investigated using the density functional theory (DFT). 
Unlike the bulk alloy, the Co-rich clusters are found to be ferromagnetic and the magnetic moment 
increases with Mn-concentration, and is larger than the moment of pure Con clusters of same size. 
For a particular sized cluster the magnetic moment increases by 2 /is/Mn-substitution, which is 
found to be independent of the size and composition. All these results are in good agreement with 
recent Stern-Gerlach (SG) experiments [Phys. Rev. B 75, 014401 (2007) and Phys. Rev. Lett. 98, 
113401 (2007)]. Likewise in bulk yinxCoi-x alloy, the local Co- moment decreases with increasing 
Mn-concentration. 

PACS numbers: 75.75.+a, 36.40.Cg, 73.22.-f 



I. INTRODUCTION 

Nanoalloy clusters have received considerable atten- 
tion for their peculiar catalytic, optical, magnetic, elec- 
tronic, and geometric properties. ^^^^^^^ For such clus- 
ters, chemical and physical properties can be tailored by 
varying not only the size but also the composition for 
a specific purpose. This opens the way to a large vari- 
ety of potential applications in areas su ch as high-density 
recording,^ catalysis j ^ ' ^'^ ' optics j-^^ l -^-^l-^^ l and biomedical. 
Particularly, the interest in transition metal clusters 
arises from a desire to seek a solution to the techno- 
logically important question: how magnetic properties 
change in reduced dimensions? The candidates chosen 
for the present study, Mn and Co have very interesting 
properties in low dimensions. Manganese, though anti- 
ferromagnetic as bulk, shows finite magnetic moment in 
reduced dimension^^^^^^ whereas. Cobalt sho ws en- 
hanced magnetic moment compared to the bulk.^^'ISEol 
Therefore, it will be interesting to see how the properties 
of the bimetallic cluster formed out of these two elements 
change with composition, atomic ordering, and size. 

The first step to the study of cluster properties is 
the determination of the ground state structures and 
the complexity to locate that increases with the clus- 
ter size as the number of local minima in the potential 
energy surface increases. This leaves a number of pos- 
sible geometric and/or magnetic isomers in pure clus- 
ters for each size. Compared to the pure clusters, in 
alloy clusters "homotops"^ are possible in addition to 
such usual isomers. Jellinek introduced the term "ho- 
motop" to describe A^^B^ alloy cluster isomers, with a 
fixed size {x + y=n) and composition {x/y ratio), which 
have the same geometrical arrangement of atoms, but 
differ in the way A- and B-type atoms are arranged. 
Due to the presence of these homotops, there arises a 



large number of combinatorial possibilities which makes 
the finding of lowest energy structures for alloy clusters 
even more computationally expensive task than that for 
pure clusters. Thus most of the theoretical studies done 
on bimetallic alloy clusters take resort to some empir- 
ical ma nybody potentia l to reduce the computational 
expense.P' ^-^ ' ^^ ' ^^l^^ -^'^^' Our study is one of the very few 
which uses an ab-initio methodology for transition metal 
nanoalloy clusters. Here we study bimetallic Mn^^Co^ 
clusters in the size ranging from 2 to 5 for all possible 
stoichiometry. Using first-principles DFT we study the 
evolution of structural, electronic and magnetic proper- 
ties as we change the size and composition. The interplay 
between these properties yield many interesting features. 
Such features are analyzed in greater depth through den- 
sity of states and partial charge density. 

In earlier works, we have studied pure Mun (Refs.E^ 
and Con (Ref.^^) clusters. A transition from ferromag- 
netic to ferrimagnetic Mn-Mn coupling is observed at n = 
5 for Mn^ clusters and the ferrimagnetic states continue 
to be the ground state for larger clusters. On the other 
hand, pure Co^ clusters are found be ferromagnetic with 
magnetic moment higher than the bulk. Calculated mag- 
netic moments of pure Mn^ and Co^ clusters show very 
good agree ment with the Stern-Gerlach molecular beam 
experiments.EHIEBIini 

Neutron scattering studies of bulk Mn^^Coi-a; alloy 
have been used to determine the variation of individ- 
ual atomic moments, /iMn and /ico, with increasing 
MniCo ratio.'^^'^^ On the Co-rich side, /iMn and /ico 
are aligned antiferromagnetically, with the magnitude of 
both /ico and /iMn decreasing monotonically with increas- 
ing Mn content such that the mean per-atom moment (/i) 
of the alloy also decreases strongly with Mn content. In- 
fact, /i decreases from 1.72 /i^ for x=0 to /i^ for x=0.32 
(Ref.^^). However, small Mn-Co clusters have been found 
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to behave in a completely different way.^^^^ Recently, 
Knickelbein has found that unlike bulk Mn^^Coi-a; al- 
loys, in which the presence of Mn forces the mean per- 
atom moment to decrease, the significant presence of Mn 
in medium sized Mn^^Co^y {n=x + ?/=ll— 29) clusters re- 
sults in overall magnetic moment that are comparable 
to those of the corresponding pure Co^ clusters, and in 
some cases (e.g. n=ll-14) even larger. More recently Yin 
et at. have measured the magnetic moments of Mn^^Co^ 
{y <60; x<y/3) clusters and found an increase in per- 
atom magnetic moment for Co-rich Mn^^Co^ cluster with 
increasing Mn concentration.^^ This enhancement in mo- 
ment due to Mn doping is independent of cluster size 
and composition. On the other hand, for Mn-rich clus- 
ters, for more than 40% Mn concentration, the average 
magnetic moment of Mn^^Co^ cluster decreases with in- 
creasing Mn concentration. This suggests that unlike 
bulk Mn-Co alloys, both Mn and Co within the small 
Mn-Co clusters retain substantial moments even at high 
Mn fractions. However, the magnitudes of individual 
Mn and Co moments could not be measured, and con- 
sequently the nature of magnetic coupling can not be 
concluded in a SG experiment. Motivated by these re- 
cent SG experiments,^^ in the present paper, we study 
Mn^^Co^y clusters from first-principles in the size range 
n=2-5. In order to make a direct comparison to the 
experiments we also study Mn2Coii cluster, which lies 
within the experimental regime 



II. COMPUTATIONAL DETAILS 

Calculations are performed usin g D FT based pseu- 
dopotential plane wave method.^^ We have cho- 
sen the projector augmented wave method, and 
used the Perdew-Bruke-Ernzerhof exchange-correlation 
functionaP^i for the spin-polarized generalized gradient 
correction. The 3d and 45 electrons of Mn and Co are 
treated as valence electrons. The wave functions are ex- 
panded in a plane wave basis set within 270 eV kinetic 
energy. Reciprocal space integrations are carried out at 
the r point. Symmetry unrestricted optimizations (of 
both geometry and spin) are performed using the con- 
jugate gradient and quasi-Newtonian methods until all 
the force components are less than a threshold value of 
0.005 eV/A. Simple cubic supercells are used with peri- 
odic boundary conditions, and it is made sure that two 
neighboring clusters are separated by at least 10 A vac- 
uum space. This ensures that the interaction of a cluster 
with its periodic image is negligible. Earlier we have 
used sa me m ethodology to study pure Mn^ and Co^ 
clusters .^^^^ For each cluster size all possible "homo- 
tops" have been considered for several geometric struc- 
tures with all possible compositions. We have also con- 
sidered all possible spin multiplicities for each of these 
structures. These ensures the robustness for the ground 
state search. It should be mentioned here that calcula- 
tions are done within the collinear spin assumption. 



The binding energy per atom (Eb) is defined as, 

EB{Mn^Coy) = - [xE{Mn) + yE{Co) - E{Mn^Coy)], 

(1) 

where x (y) is the number of Mn (Co) -atoms in 
Mn^^Co^ cluster, n(= x ^ y) is the cluster size and 
EiMiixCoy)^ £^(Mn), and ^(Co) are the total energies 
of Mn^^Co^y cluster, and an isolated Mn- and Co- atom, 
respectively. For a given n and for a certain composition, 
the structure with the highest binding energy is consid- 
ered to be the ground state. The local magnetic moment 
/ix at X-atom can be calculated from, 

/ix= / [pT(r)-Pi(r)]^r (2) 

where p| (r) and py (r) are spin- up and spin-down charge- 
densities, respectively and R is the radius of the sphere 
centered on the atom X. For a particular cluster, R is 
taken such that no two spheres overlap i.e. R is equal to 
half of the shortest bond length in that cluster. 

III. RESULTS AND DISCUSSIONS 
A. Ground states and significant isomers 

It is necessary to carry out calculations for not only 
the ground state, but also for the low energy isomers, i.e. 
clusters with different geometries, homotops and differ- 
ent magnetic arrangements, which have energies close to 
that of the ground state. This is because when, in our 
earlier works,!!^!^ theoretical results were compared with 
experimental results, it was noted that for a particular 
size of cluster, the isomers with different magnetic mo- 
ments are likely to be present in the SG beam with a 
statistical weight and essentially the measured magnetic 
moment is the weighted average of the moments of all 
those isomers. Previously we have extensively studied 
the pure Mn^ (Ref.^^) and Co^ (Ref.*^) clusters within 
the same theoretical methodology. Therefore, we do not 
elaborate pure clusters here, rather we refer the readers 
to Ref.i^ and Ref.^o. 

MuxCoy (x -\- y = 2): Due to the half-filled 3d and 
filled 4s states and due to high As^3d^ As^3d^ pro- 
motion energy Mn-atoms bind very weakly when they 
are brought together. The binding energy of Mn2 dimer 
is 0.52 eV/atom and the bond length is comparatively 
higher than all other 3d transition metal dimers.^^ On 
the other hand, C02 dimer has much higher binding en- 
ergy (1.45 eV/atom) and smaller bond length (1.96 A).!^ 
Both the pure dimers are ferromagnetic with 10 (Mn2) 
and 4: jiB (C02) moments. The bond length of MnCo 
dimer is in between the pure dimers and the binding en- 
ergy increases monotonically in Mn2<MnCo<Co2 order 
(Table The Mn-Co coupling is also ferromagnetic in 
the mixed dimer. 
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TABLE I: Structure, binding energy Eb, relative energy to 
the ground state (GS), AE=E—Egs, per-atom magnetic mo- 
ment fi and average bondlength (Lb) for the ground states of 
MuxCoy clusters with x -\- y — 2. 

Cluster Structure Eb AE /2 (Lb) 
(eV/atom) (eV) (/xs/atom) (A) 

"Co^ linear L45 OOO 2 1.96 

MnCo linear 1.09 0.00 3 2.06 

Mns linear 0.52 0.00 5 2.58 



MuxCoy (x -\- y = 3): The linear and triangular struc- 
tures were taken as initial guess and the results are shown 
in Figjl] and Table [llj The ground state structure are 
found to be triangular for all the compositions. The pure 
C03 has two degenerate ground states. The magnetic mo- 
ment is enhanced by 2 and 4 /i^ from the two degenerate 
ground state structures of C03 due to single Mn dop- 
ing. For Mn2Co the lowest energy state with a stable 
spin distribution^ is ferrimagnetic, with one of the Mn 
atoms having opposite spin to the other Mn and Co. The 
closest isomer to this is a linear ferromagnetic structure. 
However, it lies 0.39 eV above the ground state. Here the 
Mn-Mn distance is large (4.54 A) causing them to be fer- 
romagnetically coupled. The ferromagnetic ground state 
of pure Mna was found to be nearly degenerate with a 
ferrimagnetic structure.'!^ 




(3, 0.00 eV, 5|Ib) (3, 0.00 eV, 




(2, 0.00 eV, 9|Ib) (1, 0.00 eV, Slig) (1, 0.13 eV, III^b) 




(0, 0.00 eV, ISlig) (0, 0.05 eV, 5[i^) 

FIG. 1: (Color online) Ground state structures and significant 
isomers of Mn^cCo^ (x + y = 3) clusters. Lighter (Magenta) 
shades represent Co-atoms and darker (blue) shades represent 
Mn- atoms. The bond lengths are quoted in A. The numbers 
in the parenthesis represent the number of Co-atoms (y), dif- 
ference in total energy from the corresponding ground state 
(AE), and total magnetic moment, respectively. Down arrow 
represents the atom with antiparallel moment. 



TABLE II: Same as Table |T] for x + y^3 clusters. 



Cluster 


Structure 


Eb 












(eV/atom) 


(eV) (m 


B 1 atom) 


(A) 


C03 


Triangular 


1.78 


0.00 


1.67 


2.16 




Triangular 


1.78 


0.00 


2.33 


2.19 


MnCo2 


Triangular 


1.67 


0.00 


3.00 


2.26 


Mn2Co 


Triangular 


1.16 


0.00 


1.00 


2.40 




Linear 


1.12 


0.13 


3.67 


2.27 




(with endon Mn) 










Mns 


Triangular 


0.82 


0.00 


5.00 


2.74 




Triangular 


0.80 


0.05 


1.67 


2.48 



MuxCoy (x -\- y = 4): For these clusters square and 
tetrahedral structres were considered. Results are sum- 
marized in Table IIIII and some of the relaxed structures 
are shown in Figj2j The ground state structures come 
out to be tetrahedral for all compositions. C04 has a 
Jahn- Teller distorted tetrahedr al g round stat^ similar 
to that obtained by Castro et a/.p^and the closest isomer 
is a rhombus. Single Mn-substitution (MnCos) increases 
the ground state moment by 2 jub- The first isomer is a 
planar structure, a quadrilateral, having equal magnetic 
moment as the tetrahedral ground state. The next two 
closest isomers also have tetrahedral structure: one with 
average magnetic moment 0.5 /i^/atom lower and aver- 
age bond length 0.03 A shorter and the other with aver- 
age bond length 0.01 A longer and average magnetic mo- 
ment 0.5 /i^/atom higher than that of the ground state. 
Thus the isomer with higher bond length has higher mo- 
ment. This trend is seen for all compositions and sizes 
(Tables I to IV). The closest non- ferromagnetic MnCoa 
isomer lies much higher (0.82 eV) than the ground state 
and has zero net moment. 

Substitution of another Co with Mn (Mn2Co2) further 
enhances the total magnetic moment by 2 /i^ compared 
to MnCoa cluster. In this case, the closest ferrimagnetic 
isomer, with average magnetic moment 1 /i^/atom, lies 



TABLE III: Same as Table [I] for X -\- y — 4 clusters. 



Cluster 


Structure 


Eb AE 
(eV/atom) (eV) (/i^ 


; / atom) 


(Lb) 

(A) 


C04 


Tetrahedral 


2.28 


0.00 


2.5 


2.34 




Rhombus 


2.25 


0.11 


2.5 


2.25 


MnCos 


Tetrahedral 


2.06 


0.00 


3.0 


2.34 




Quadrilateral 


2.05 


0.05 


3.0 


2.32 




Tetrahedral 


2.02 


0.18 


2.5 


2.31 




Tetrahedral 


2.00 


0.25 


3.5 


2.35 


Mn2Co2 


Tetrahedral 


1.87 


0.00 


3.5 


2.41 




Tetrahedral 


1.82 


0.22 


4.0 


2.46 


Mns Co 


Tetrahedral 


1.60 


0.00 


4.0 


2.54 




Tetrahedral 


1.54 


0.22 


4.5 


2.60 




Rhombus 


1.51 


0.35 


1.5 


2.39 


Mn4 


Tetrahedral 


1.18 


0.00 


5.0 


2.54 




Tetrahedral 


1.16 


0.08 


2.5 


2.59 




Tetrahedral 


1.13 


0.20 


2.0 


2.58 




Tetrahedral 


1.13 


0.20 


0.0 


2.65 



4 




(4, 0.00 eV, 10 lag) (4, 0.11 eV, 10 [i^) (3, 0.00 eV, 12 




(3, 0.05 eV, 12 (2, 0.00 eV, 14 ^ig) (1, 0.00 eV, 16 [i^) 




(1, 0.35 eV, 06 jig) (0, 0.00 eV, 20 jig) (0, 0.08 eV, 10 iig) 

FIG. 2: (Color online) Ground state structures and significant 
isomers of Mn^cCoy {x -\- y = 4) nanoalloy clusters. Same 
conventions are followed as in Fig[l] 

far above (0.56 eV) from the ground state. A further 
substitution of a Co- atom with Mn atom (MnaCo) again 
leads to further 2 /i^ increase in magnetic moment. A 
planar (rhombus) ferrimagnetic isomer lies 0.35 eV above 
the ground state. Mn4 has binding energy 1.07 eV/atom, 
which is much lower than that of C04 and is ferromag- 
netic. However, it has a close ferrimagnetic isomer, which 
lies only 0.08 eV higher in energy. It is interesting to note 
here that for these clusters with n=4, the energy split- 
ting between the ferromagnetic ground state and optimal 
ferrimagnetic state decreases with increasing Mn. 

MuxCoy (x -\- y = 5): The trigonal bi- pyramidal 
(TBP), square pyramidal (SQPD) and (planar) pentago- 
nal geometries are considered as starting configurations. 
The results are tabulated in Table IIVI and minimum en- 
ergy structures are shown in Fig|3] All the clusters with 
different compositions have trigonal bi-pyramidal ground 
state. For MnCo4 cluster, the ground state magnetic 
moment is enhanced by 2 /i^ compared to the pure Co- 
pentamer due to single Mn substitution. The single-Mn 
doping also enhances the energy splitting between TBP 
and SQPD structures. We have found a similar trend in 
magnetic moment for Mn2Co3 as it is further enhanced 
hy 2 jiB due to another Mn-substitution. For this cluster 
the nearest 'homotop' lies 0.15 eV higher. 

The incremental behavior of magnetic moment with 
Mn-substitution is no longer observed as the clusters be- 
come ferrimagnetic, when the number of Mn- atoms in- 
crease to three or more, i.e., when the clusters become 
Mn-rich. Thus, magnetic moment of Mn3Co2 drops by 4 
jiB from the pure C05, as in this cluster one of the Mn- 



TABLE IV: Same as Table [T] for x + 2/ = 5 clusters. 



Cluster 


Structure 


Eb 






{Lb) 






(eV/atom) (eV) (/is 


; / atom) 


(A) 


C05 


TBP 


z.oo 


U.UU 


2.6 


2.34 




SQPD 


Z.OZ 


U.lo 


2.6 


2.27 


MnCo4 


TBP 


Z.4D 


n nn 

U.UU 


3 


2.38 




SQPD 


9/11 


n 9Q 
U.zo 


2.6 


2.29 




SQPD 


9 /in 
Z.4U 


n 9Q 
U.zy 


3 


2.35 


Mn2Co3 


TBP 


9 91^ 
Z.ZO 


n nn 
U.UU 


3.4 


2.43 




TBP 


9 9/1 
Z.Z4 


n nQ 
U.Uo 


3 


2.40 




TBP 


9 99 

z.zz 


n 1 
U.lo 


3.4 


2.42 




TBP 


9 99 
Z.ZZ 


n 1 
U.lo 


3.4 


2.37 


Mn3Co2 


TBP 


z.uo 


D nn 

U.UU 


1.8 


2.45 




TBP 


9 ni 

Z.Ui 


n 1 n 
U.lU 


1.8 


2.45 




TBP 


1 QQ 


n 1 Q 
u.iy 


3.8 


2.50 




TBP 


1 QQ 


n 99 
u.zz 


3.4 


2.44 




TBP 


1 Q7 


n Qn 
U.oU 


3.8 


2.54 


Mn4Co 


TBP 




n nn 
U.UU 


2.2 


2.52 




TBP 


1.76 


0.02 


0.2 


2.49 




TBP 


1.75 


0.04 


4.2 


2.55 




TGPD 


1.73 


0.15 


2.2 


2.52 




SQPD 


1.72 


0.18 


3.8 


2.51 




TBP 


1.71 


0.25 


4.2 


2.60 




TBP 


1.70 


0.30 


1.8 


2.51 


Mns 


TBP 


1.41 


0.00 


0.6 


2.54 




TBP 


1.40 


0.06 


2.6 


2.58 




TBP 


1.40 


0.07 


1.0 


2.59 



atoms has antiparallel spin alignment with the others. 
The first isomer found is a homotop of the ground state. 
A ferromagnetic homotop, which follows an increment of 
2 /i^ /Mn-substitution, lies much higher (0.30 eV), and 
has larger (~ 2.80 A) Mn|-Mn| separations. In contrast, 
in the ferrimagnetic ground state Mn|-Mn| distance is 
much shorter (2.44 A) than the Mn|-Mn| distance (2.88 
A). This is in general true for other clusters (Figj3|. The 
Mn4Co is also ferrimagnetic, which has 2/1^ less moment 
than that of pure C05. However, it has a ferromagnetic 
structure which obeys the '2 /i^ /Mn-substitution incre- 
ment' rule and this structure lies only 0.04 eV higher. 
This cluster also has two different homotops in the form 
of a distorted tetragonal pyramid (TGPD) and a SQPD 
and they lie 0.15 and 0.18 eV higher, respectively. It has 
previously been reported that pure Mns is ferrimagnetic 
in its ground state and a ferromagnetic isomer lies 0.19 
eV higher.-^ 



B. Binding energy and stability 

The coodination number increases with cluster size, 
and thus the binding energy, which we have discussed 
earlier in detail for pure Mn^ and Co^ clusters. iJXM 
However, the binding energy of pure Co^ clusters are 
much larger than that of pure Mn^ clusters of same size. 
This is because the Mn-atoms have half-filled 3(i and filled 
As shells and also have high 4s^3(i^ As^'id^ promotion 
energy. As a consequence, the Mn-atoms do not bind 
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(5, 0.00 eV, 13 |LiB ) 




2.74 

^ ^.41 

(3, 0.00 eV, 17 ILiB ) 




2.11 



(2, O.lOeV, 9 |LIb) 

2.53 




2.73 

(l,0.15eV, 11 ILIb) 




2.17 

(5,0.15eV, 13 |LIb) 





2.72 

(l,0.18eV, 19 ILIb) 




(4, 0.00 eV, 15 ILIb ) 





(l,0.00eV, 11 ILIb) 



2.68 



(1,0.25 eV, 21 ILIb) 




2.81 




2.21 

(4, 0.23 eV, 13 ILIb) 




(3, 0.03 eV, 15 ILIb ) (3, 0.15 eV, 17 |LIb ) (3, 0.15 eV, 17 ILIb ) 




2.57 



2.80 



(l,0.02eV, 1 |LIb) 

\^ 2.85 




(l,0.30eV, 9 I^b) 




(0, 0.06 eV, 13 ILIb ) 



(0,0.07eV, 5 ILIb) 



2.15 

(4, 0.29 eV, 15 ILIb ) 



2.70 




(0, O.OOeV, 3 |LIb) 



FIG. 3: (color online) Ground state structures and significant isomers of Mn^Coy (x + ^ = 5) nanoalloy clusters. Same 
conventions are followed as in Fig^ 



strongly when they are brought tog 

ether™ 

The Figj4] shows that, in general, the binding energy 
decreases with increasing number of Mn-atoms (x) for 
a particular sized cluster n. This again is due to the 
weaker binding among Mn atoms, and relatively weak 
Mn-Co binding than Co-Co bonding. 

The stability S of these nanoalloy clusters with the 
variation of Mn-atoms x in it can be defined as, 

S{n, x) = E{n, x + 1) + E{n, x - 1) - 2E{n, x), (3) 

where E{n^x) is the bound state energy of the Mn^^Co^ 
{n = x^y) cluster, and is shown in Fig js] for n = 5. The 
MnCo4 clusters is seen to have maximum stability. A 
sudden dip in stability is seen for Mn2Co3, which can be 
described in terms of different bond distributions. We see 




FIG. 4: (color online) Binding energy of MuxCoy alloy clus- 
ters as a function of Mn-atoms. 
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2.15 




2.05 



1.95 



X in Mn^COy 



X in Mnjj.COy 



12 3 
X in Mnjj.Co2 



I 

o 




12 3 
X in Mnjj.Co2 



FIG. 5: (color online) (a) Relative stability and (b) different 
(Mn-Mn, Mn-Co, and Co-Co) bond lengths of n = 5 cluster 
as a function of Mn- atoms. 

that Mn2Co3 has maximum (minimum) Mn-Mn (Co-Co) 
bond length, and thus, has the weakest (strongest) Mn- 
Mn (Co-Co) interaction among all the n=5 clusters with 
different compositions. This reduced hybridization of the 
Mn-atoms in turn reduces the stability of this cluster 
compared to its neighbors in x, and also leaves the clus- 
ter ferromagnetic. Similarly, for n=4 cluster, Mn2Co2 
has larger average Mn-Mn distance (2.83 A) than MnsCo 
(2.74 A) and found to be less stable. 

C. Magnetic moment 

The atoms in the pure Co^ clusters are found to be 
aligned ferromagnetically,^^ whereas a ferromagnetic to 
ferrimagnetic transition has been observed for pure Mn^ 




' ' ' ' ' ^ 

12 3 4 

X in Mn_j.Coy 



FIG. 6: (color online) Calculated magnetic moment per atom 
as a function of Mn-atoms. 



FIG. 7: (a) Co-Co bond length and (b) average local Co- 
moment (/ico) in MnajCo2 clusters as a function of x. 



at n = 5 and remains the same for larger sized clusters .^^^ 
The local magnetic moment of Mn {fiun) is found to 
be higher than /ico due to more number of unpaired d- 
electrons in Mn-atom {Sd^As'^) than in Co- atom {Sd^As'^). 

Fig|6] shows the average magnetic moment /i as a func- 
tion of the number of Mn-atoms, x, for different cluster 
sizes n. For n = 2 and 4, ft increases with increasing 
number of Mn-atoms. This is expected since for these 
sizes, all the clusters with different compositions are fer- 
romagnetic. As mentioned earlier, the net number of un- 
paired d-electrons in Mn (five) is greater than that in Co 
(three). Thus as Co-atoms are replaced by Mn-atoms, to- 
tal cluster moment increases by ~ 2 /i^/Mn-substitution. 
For cluster size n=3 and composition x=2, and also for 
n=5 and x > 3, i.e., as the clusters become Mn-rich, 
ft decreases with x. We see that these Mn-rich clusters 
are ferrimagnetic. For pure Mn^ clusters ferrimagnetic 
coupling is observed for n > 5 (Ref.^^). On the other 
hand, for bimetallic Mn^^Co^y clusters, ferrimagnetic cou- 
pling of Mn-atoms is observed for n=5 with x > 3. This 
could be because, for pure Mns and Mn4, the ferrimag- 
netic isomers having average magnetic moment 1.67 and 
2.50 /i^/atom, respectively, lie just 0.05 eV and 0.08 eV 
above their corresponding ferromagnetic ground states. 
Thus perturbing these clusters with Co-atoms could in- 
duce the corresponding ground state to be ferrimagnetic. 

Another interesting observation is made that if we dope 
a C02 dimer with increasing number of Mn-atoms, the 
Co-Co bond length k.pdf increasing (FigjT^). This indi- 
cates that Mn-atoms reduce the hybridization between 
Co-atoms in the dimer. Therefore, we have further stud- 
ied the average /ico iii C02 dimer with increasing x, 
and is shown in Figj7]3. It is seen that /ico remain al- 
most the same for all the Mna^Co2 clusters (1.94, 1.88, 
1.88, and 1.85 /i^/atom for C02, MnCo2, Mn2Co2, and 
Mn3Co2, respectively). Although the Co-Co bond length 
in Mna^Co2 increases with increasing x, the coordination 
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FIG. 8: (Color online) Isosurfaces of 
magnetization density for (a) Mn2Co, 
(b) MnaCo, and (c) Mn4Co corre- 
sponding to 0.02, 0.03, and 0.04 e/A^ 
respectively. Dark (Blue) and light 
(magenta) color represent Mn- and 
Co- atoms, respectively. The Mn^ is 
marked by down arrow. Blue (Red) 
surface indicates positive (negative) 
magnetization density. 



of Co-atoms also increases as we go along C02 MnCo2 
Mn2Co2 Mn3Co2, which leaves the effective hy- 
bridization of the Co-atoms unaffected. 

The magnetization density (difference between the up 
and down spin densities) further illustrates (Figj8| the 
magnetic nature of these nanoalloy clusters. It is interest- 
ing to note that in ferrimagnetic Mn2Co (Mn4Co) clus- 
ter the Mn|-Mn| separation is 2.64 A (2.61 and 2.65 A). 
This is significantly smaller than Mn^-Mn^ separation 
in ferromagnetic Mns and Mn4 clusters, 2.74 and 2.70, 
respectively.^^ Also the average Mn|-Mn| bond in ferri- 
magnetic Mn4Co is longer (Figjsj^) than the Mn|-Mn| 
one. Thus it seems that wherever there is a contraction 
of Mn-Mn bond due to Co doping, the two Mn atoms 
concerned get antiferromagnetically aligned. This is fur- 
ther supported by the fact that for MnaCo, where the 
distances between the Mn atoms (2.74 A) is not reduced 
by Co-doping, the Mn atoms remain ferromagnetically 
coupled. The dependence of Mn-Mn coupling on the sep- 
aration is due to the modification of hybridization. We 
will further discuss the magnetic nature of the clusters in 



Section IIIH through partial charge density. 



sizes where an exact 50-50 concentration is not possible. 
We see a 80-31% enhancement in magnetic moment com- 
pared to pure Con clusters. As we have already seen that 
these nanoalloy clusters are ferromagnetic unless the clus- 
ter is Mn-rich, thus when we replace Co-atoms in pure 
Co^ clusters with Mn, the magnetic moment increases. 

In a separate SG experiment, Yin et al. has found 
that magnetic moment increases as the number of Co- 
atoms, in Mn^^Co^y increases in the range 15-45 keep- 
ing X constant at 1 — 9 (i.e., all the cluster are Co-rich). 
This observation also holds good in the present study for 
smaller clusters as long as the alloy cluster is Co-rich. 
For an example, as we go along MnCo (6 /i^) MnCo2 
(9 11b) MnCo3 (12 jis) ^ MnCo4 (15 113) the to- 
tal moment increases with increasing y. They have also 
found that the total enhancement in the magnetic mo- 
ment to be 1.7 /iB /Mn-substitution, and is independent 
of the cluster size n and composition, which is justified by 
the virtual bond state model. The present calculation 
is in good agreement and we find this enhancement to be 
2 /i^ /Mn-substitution (Tables IpV) as long as cluster is 
ferromagnetic (Co-rich). 



D. Comparison of magnetic moment with 
Stern- Gerlach experiment 



In a recent SG experiment, Mn^^Co^ (n=ll-29) clusters 
were produced via pulsed laser vaporization from cylin- 
drical MiixCoi-x alloy (Mno.15Coo.85 or Mno.5Coo.5) tar- 
get rods.^^ The mass spectrum shows that average Mn- 
fraction in the alloy clusters are always less than in the 
corresponding source rod. However, this Mn-fraction in- 
creases with the increase in cluster size. Therefore, in the 
experiment, alloy clusters were always Co-rich compared 
to the bulk source. Clusters generated from both the 
Mno.15Coo.85 and Mno.5Coo.5 target rods have moments 
that are similar to that of the pure cobalt clusters of same 
size. In fact in smaller size range, n=ll-14, the average 
moment is enhanced by 28-60% for clusters generated 
from Mno.5oCoo.5o source. We have compared the mag- 
netic moments of pure Co^ clusters with 50-50 Mn-Co 
nanoalloy clusters (Figj9| in the size range n =2-5 to see 
whether the observed experimental trend continues for 
the smaller sizes. It should be pointed out here that, 
for Mno.5oCoo.5o clusters, we have considered the Co- 
rich clusters to mimic the experimental situation, for 
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FIG. 9: (color online) Magnetic moment of pure Con and 
Mno.sCoo.s clusters. For cluster sizes n = 3 and 5, the cobalt 
rich clusters (MnCo2 and Mn2Co3, respectively) are consid- 
ered. See the text. 
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FIG. 10: (color online) Average local Co- moment (/xco) as a 
function of Mn-atoms for different cluster sizes. 



E. Comparison with bulk alloy 

The individual atomic moments of Mn and Co for 
Mna;Coi_a; bulk alloy have been determined with increas- 
ing MniCo ratio through neutron scattering studies^^. 
On the Co-rich side /iMn and /ico are aligned antiferro- 
magnetically, with the magnitude of both /iMn and /ico 
decreasing monotonically with increasing Mn-content. 
Same trend is also seen for the average magnetic moment 
of the alloy fl=xjiiMn +(1 — x) /ico- However, it is seen 
from the SG experiment d^^ * ^^ * that Co-rich Mn^^Coi-a^ 
nanoalloy clusters retain substantial magnetic moment 
in the size range n= 11-29. Moreover, for clusters with 
n=ll-14 obtained by laser vaporization of Mno.15Coo.85 
rods, the magnetic moment increases by 88-148% com- 
pared to the bulk value. At 20% Mn-concentration for 
cluster size n=5 we estimate 241% enhancement in mag- 
netic moment from the corresponding bulk value. This 
enhancement of p in low dimension over bulk alloy is 
not only due to decrease in coordination number, which 
effectively reduces the hybridization among the orbit als, 
but also due to the ferromagnetic Mn-Co coupling (unlike 
bulk alloy) in Co-rich Mn^^Co^y nanoalloy clusters. 

The average fiun and fico (calculated using Eq. 2) 
are studied individually for different sizes with increas- 
ing Mn-concentration. The average /iMn oscillates with 
increasing x. For clusters like Mn2Co which are ferrimag- 
netic, the average magnetic moment of Mn is very low 
(0.08/iB/atom), whereas for ferromagnetic (like Mn2Co2) 
clusters it is high (4.21/iB/atom). Here, it is to be noted 
that even for such ferrimagnetic clusters the individual 
/iMn values are high but due do their ferrimagnetic Mn- 
Mn alignment the net Mn-moment is low. The average 
/iCo is plotted in FigjToj It is clear from Fig 10 that. 




likewise in bulk alloy, the average /ico decreases with 
increasing Mn concentration for all cluster sizes,'^ and, 
as we will see in Section HI G[ is because of increasing 
Mn neighbor. However, /iMn doesn't behave in similar 



(35 |^B,0.15eV) 

FIG. 11: (Color online) The ground state (a), first (b), and 
second (c) isomers for CoiiMn2 nanoalloy cluster. Same color 
conventions are used as in Figjl] The numbers in the paren- 
thesis represent the total moment and relative energy to the 
ground state, A^, respectively. 



fashion, and unlike the bulk alloy, in Co-rich Mn^^Co^ 
nanoalloy clusters, /iMn and /ico are ferromagnetically 
coupled. 



F. A medium sized cluster (n =13) and direct 
comparison with SG experiment 

In order to perform a direct comparison with the 
SG experimeniP3 ^e have studied a Mn^^Co^ cluster of 
size n=13, with 15% Mn-concentration i.e., Mn2Coii 
nanoalloy cluster. The ground state is found to be 
nearly degenerate with magnetic moments 1.92 and 2.38 
/i^/atom. The ground state has binding energy 3.18 
eV/atom (FigjTT^), whereas the nearly degenerate isomer 
lies only 0.01 eV higher (Figjllb). In the first structure 
the Mn-Mn distance is 2.46 Aand one of the Mn-atom 
sits at the center, and consequently the magnetic mo- 
ment of this (Mn) atom gets quenched. Whereas, the 
Mn-Mn distance is comparatively larger (4.78 A) in the 
second cluster as both the Mn-atoms sit on the surface. 
These are the reasons that the second structure has larger 
magnetic moment than the first one. Another icosahe- 
dral isomer (FigjTl]^) has large magnetic moment of 2.69 
/i^/atom and lies 0.15 eV higher than the ground state. 

The icosahedral ground state structure for this Mn- 
doped cluster is quite interesting. For pure C013 the 
stable ground state is a distorted hexagonal structur^^ 
and, on the other hand, we have previously found an 
icosahedral ground state for pure Mni3 cluster.!!^ Re- 
placing two Co-atoms in C013 with Mn results in a 
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FIG. 12: (Color online) The d-projected density of states for 
the (a) central, and (b) surface Mn-atom of Mn2Coii clus- 
ter in its ground state (marked as 1, and 2, respectively, in 
Fig|ll^). A Gaussian broadening of 0.1 eV has been used. 
Corresponding integrated density of states up to the Fermi 
level, Ef, are also shown for the (c) central, and (d) surface 
Mn-atoms. 



structural change from a distorted-hexagonal-like struc- 
ture to an icosahedron. The calculated magnetic mo- 
ment of the pure C013 is found to be 25 /iB, which is 
in good agreement with experimental values, 26-30 jj^b 
(Refs.*^^^^. On the other hand, two degenerate ground 
states of Mn2Coii nanoalloy cluster have 25-31 /j^b mag- 
netic moments, which is 0-6 /ib larger than that of (cal- 
culated) the pure C013 cluster. This is in good agree- 
ment with the experiment: For n=13 cluster produced 
from Mno.15Coo.85 bulk source an enhancement in aver- 
age magnetic moment is seen over the corresponding pure 
C013 cluster. It should also be pointed out here that the 
SG experiment was done at finite temperature (91 K), so 
the cluster beam may contain both the degenerate states 
which are separated only by 0.01 eV energy. 



G. Projected density of states 

To understand the decrease in average /ico when the 
environment is made more Mn-rich than Co (Fig{To|, the 
integrated projected density of states (IPDOS) of the 
d-orbital has been investigated for the cobalt atoms in 
Mn2Co3 and Mn4Co clusters. For the first cluster in the 
TBP ground state (Figj3|, the IPDOS of the Co-atom, 
which is co-planer with two Mn-atoms, is 4.41 for the 
majority spin channel, whereas it is 2.30 for the minor- 
ity spin. For a off-plane Co-atom in same structure, the 
IPDOS is 4.30 for the majority spin channel and 2.40 
for the minority spin one. The reason for this greater 
spin polarization of the co-planar Co-atom is twofold. 



Firstly, the in-plane Co-atom has more (two) Co-atoms 
as its nearest neighbors than the off-planar ones (one). 
Secondly, the in-plane Mn-Co distance (2.62 A) is larger 
than that for the off-planar Mn-Co distance (2.41 A). 
Thus the effect of Mn-atoms felt by the in-plane Co-atom 
is less than that by the off-planar ones. This reinforces 
the observation that a Co-rich environment favors more 
spin polarization in Co than a Mn-rich environment. For 
the Co-atom in the second cluster, Mn4Co, the IPDOS 
for the majority (minority) spin channel is 3.81 (2.50), 
giving rise to a decrease in spin polarization (compared 
to average d-polarization of Co- atoms in Mn2Co3), when 
the environment is made more Mn-rich at the cost of Co. 

The d-projected density of states (d-DOS) and the cor- 
responding IPDOS for Mn2Coii are also studied. These 



are shown in Fig 12 for the central, and apex Mn-atoms 
labeled as 1 and 2, respectively, in FigfTT^. We see that 
for the central Mn-atom (Fig 12 i), the spin distribution 



in the majority and minority spin channels are nearly 
equal. Consequently, the magnitude of IPDOS (Fig|T2|3) 
for the majority and minority spin channels are close, 
2.75 and 2.35, respectively. Thus for this central Mn- 
atom spin polarization is very low. On the other hand for 
the surface Mn-atom (Fig|T2]i), the majority spin chan- 
nel is almost fully occupiecr(Fig 12 3) and the magnitude 



of corresponding IPDOS for the majority and minority 
spin channels are 4.32 and 0.61, respectively. This gives 
a large spin polarization on this atom. This is because 
the surface Mn is less coordinated and hence has less hy- 
bridization than the central one. Further, the Co-atoms 
marked as 3 and 4 in FigfTT^ have equal coordination 
number and the former has slightly larger (2.43 A) av- 
erage bond length than the later one (2.38 A). Although 
we find that the Co-atom marked as 3 has slightly weaker 
spin polarization (1.74 /i^) than that of marked as 4 (1.83 
/hb)^ as the former is bonded with one more Mn-atom 
than the later. This once again confirms our observation 
that /ico decreases with increasing Mn neighbors. 



H. Partial charge densities 

In order to take a deeper look into the electronic be- 
havior responsible for antiparallel alignment of one fiun 
to other Co- and Mn-atoms in Mn4Co, we have inves- 
tigated the eigenvalue spectrum of the occupied levels. 
We find among the eight highest occupied levels (i.e., 
HOMO — levels, where = — 7), four of them are 
occupied by minority spins. On the other hand, among 
the eleven deepest occupied levels, six levels are occupied 
by minority electrons. Comparing this spin distribution 
picture with that of Mn3Co, we see that for this ferro- 
magnetic cluster also, only four of the eight highest occu- 
pied levels are occupied by minority spins likewise in the 
ferrimagnetic Mn4Co. However, only one of the deepest 
eleven states are occupied with minority spin. There- 
fore, it seems that these deep valence levels play a major 
role in determining the magnetic structure of these two 
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FIG. 13: (Color online) The charge density isosurfaces of the 
six down electrons in the deepest energy levels for Mn4Co 
cluster. They are shown at (a) 0.04, (b) 0.01, (c) 0.01, (d) 
0.01, (e) 0.01, and (f) 0.03 e/A^ respectively. Dark (Blue) 
and light (magenta) color represent Mn- and Co- atoms, re- 
spectively, and the Mn^ is marked by down arrow. 



(a) 





between this Mn^ and the Co-atom for lower isosurface 
values (Fig, 13). The shape of the charge densities in- 
dicates that these electrons are d-electrons. Thus, single 
Co-doping in ferromagnetic Mn4 makes Mn4Co ferrimag- 
netic. On the other hand, the charge densities of the ma- 
jority electrons are spread out among the Co- atom and 
the three Mn^ -atoms at small isosurface values. Charge 
densities for two of these majority electrons are shown in 
FiglMl The shape of these charge density isosurfaces also 



exhibits d-orbital character. It is interesting to observe 
that the Mn^ is completely devoid of majority charge 
contribution for these deep valence levels. Thus it must 
be so that the Mn^ goes from the As^Zd^ to As^Zcfi con- 
figuration in Mn4Co. Due to Co doping the Mn|-Mn| 
separation in the cluster is reduced, which consequently 
enhances the hybridization and cause the six d-electrons 
in this Mn| to belong to the minority spin channel. 

In contrast the charge density of the electrons occupy- 
ing the deepest energy levels of MnsCo cluster exhibits 
quite a different picture. The charge densities of both 
the majority and minority electrons for small isosurface 
value (Fig 15 ) are spread out amongst all the four atoms 
in the cluster and not localized on any one of them. 




t 

FIG. 14: (Color online) The charge density distribution for 
two of the majority electrons in the deep levels for Mn4Co 
cluster. Dark (Blue) and light (magenta) color represent Mn- 
and Co-atoms, respectively, and the Mn^ is marked by down 
arrow. The isosurfaces are drawn at 0.04 e/A^ density. 




FIG. 15: (Color online) The charge density distribution for 
(a) the minority electron at 0.01 e/ isodensity and (b) for 
one of the majority electrons at 0.03 e/A^ isodensity in the 
deep levels of MnsCo cluster. Dark (Blue) and light (ma- 
genta) color represent Mn- and Co- atoms, respectively. 



IV. SUMMARY 

We have studied structure, bonding, and magnetism in 
small bimetallic Mn^^Co^y (x + ?/=2-5) clusters from first- 
principles DFT calculation. Due to van der Waals like 
weaker bonding among Mn-atoms and relatively strong 
Co-Co bonding than Mn-Co bonding, the binding ener- 
gies of the alloy clusters decrease with increasing Mn- 
concentration. Interesting effects in binding energy, sta- 
bility and magnetism in the nanoalloy clusters are ex- 
plained through the interplay between bond length and 
coordination. The Co-rich clusters are found to be fer- 
romagnetic unlike the bulk alloy and the corresponding 
magnetic moment is higher than the pure C o^ clusters as 
is seen in the recent SG experiments j ^^ * ^^ ^ Moreover, the 
magnetic moment of Co-rich nanoalloy clusters increase 
with Mn-concentration and this increment is 2 /i^/Mn- 
substitution and is independent of cluster size and com- 
position. Co-atoms are found to be more magnetically 
polarized in a Co-rich environment than in Mn-rich one, 
i.e., likewise in bulk alloy, as the environment is made 
more Mn-rich, the average /ico decreases. 
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